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Alpha-particle (quartet) condensation in homogeneous spin-isospin symmetric nuclear matter is 
investigated. The usual Thouless criterion for the critical temperature is extended to the quartet 
case. The in-medium four-body problem is strongly simplified by the use of a momentum projected 
mean field ansatz for the quartet. The self-consistent single particle wave functions are shown and 
discussed for various values of the density at the critical temperature. 

PACS numbers: 03.75.Nt, 03.75.Ss, 21.65.-f, 74.20.Fg 
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I. INTRODUCTION 

Investigation of pairing in different Fermi systems is 
still on the forefront of active research. Examples are nu- 
clear physics [1| and the physics of cold fermionic atoms 
Q ■ The formation and condensation of heavier clusters 
in Fermi systems is much less studied. 

In cold atom physics the recent advent of trapping 
three different species of fermions Q has opened up the 
possibility of creating gases of heavier clusters. For the 
time being those may be trions (bound state of three dif- 
ferent fermions) but in the future one also can think of 
quartets (bound state of four different fermions). The 
latter are specially interesting because of their bosonic 
nature and the possibility of Bose-Einstein condensation 
(BEC) of quartets. The description of quartet conden- 
sation to occur has been attempted with an extension of 
the so-called Cooper problem to the four body case in 
[j). In § a variational procedure for the condensation of 
(2s + l)-component fermion clusters, with s the fermion 
spin, has been proposed. A quartet phase has been found 
in a one dimensional model with four different fermions 

On the other hand, nuclear physics, because it is 
a four-component fermion-system (proton/neutron spin 
up/down), all fermions attracting one another, leading 
to the very strongly bound a-particle, is a proto-type 
system for quartetting. There, the formation of clusters 
has been an object of study almost since the beginning 
of nuclear physics Q. Of course, pairing also exists in 
nuclei from where it is concluded that neutron stars are 
superfluid. Nuclei are very small quantum objects with 
only a (slowly) fluctuating phase (the conjugate variable 
to particle number N). Actually the number of Cooper 
pairs in nuclei generally does not exceed about a dozen 
(often much less) and yet clear signs of superfluidity are 
observed in nuclei (e.g., moments of inertia strongly re- 



duced from their classical value), implying that no critical 
size exists from where signatures of superfluidity abruptly 
disappear. One, thus, can safely extrapolate from finite 
nuclei to superfluidity in neutron stars. On the other 
hand, as already mentioned, in nuclear physics the ex- 
istence of quartets (a-particles) as subclusters of nuclei 
is omnipresent. As well known, many lighter nuclei with 
equal proton and neutron numbers [Z = N) show, for in- 
stance in excited states, strong a clustering. The concept 
that these a-particles may form a condensate in certain 
low density states of nuclei and that this may, in analogy 
to the pairing case, be a precursor sign of a-particle con- 
densation in infinite matter 0] , has come up only recently 
|9| . Also heavy nuclei seem to have preformed a-clusters 
in the surface because of their well known spontaneous a 
decay properties. 

Symmetric nuclear matter does not exist in nature be- 
cause of the too strong Coulomb repulsion. However, in 
collapsing stars, so-called proto- neutron stars, the frac- 
tion of protons is still high [lOj and the formation of a- 
particles and, at sufficiently low temperature, their con- 
densation may eventually be possible. At any rate, it 
seems evident that nuclear matter at various degrees of 
asymmetry is unstable with respect to cluster formation 
in the low density regime. At zero temperature, the most 
stable nucleus is 56 Fe but as a function of temperature, 
density, asymmetry, other cluster compositions of infi- 
nite baryonic matter may be formed. Several theoretical 
studies predict a-phases to exist in certain temperature- 
density- asymmetry domains [Til ]. 

In view of the complexity of the task, the objective 
of the present work is quite modest. We want to study 
the critical temperature of a-particle condensation as a 
function of density and temperature in symmetric nuclear 
matter. Still, even this task will not be carried out down 
to the BEC limit. We will study the critical temperature 
T" for the onset of formation of a-particles in a thermal 
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FIG. 1: Sketch of a-particle configuration, indicating that the 
two protons and two neutrons occupy the lowest OS level in 
the mean field potential of harmonic oscillator shape. 



gas of nucleons. This shall be done with a theory anal- 
ogous to the famous Thouless criterion for the onset of 
formation of Cooper pairs in a superconductor. On the 
microscopic level the problem is still very challenging, 
since it amounts to solve an in-medium four-body prob- 
lem. In spite of that, solutions have already been worked 
out in the past, either solving the Faddeev-Yakubovsky 
equations [3] or with an approximate procedure Q. 

In this work, we will continue along those lines. The 
final objective is to reach the BEC regime in a treatment 
similar to the one of Nozieres and Schmitt-Rink (NSR) 
theory [13j, but for quartets. Needless to say that this 
only will be possible if the whole formalism can radically 
be simplified. Actually, as we will show in this work, 
such a procedure may well exist. In any case, it is not 
conceivable that one treats condensation of bosonic clus- 
ters built out of N fermions on the level of non-linear 
in medium TV-body equations for N > 2. On the other 
hand, it is well known, that nuclei can satisfactorily be 
described, in mean field approximation [T^ |. Projecting 
these mean field (Hartree-Fock) type of solutions on zero 
total momentum (K = 0) will then allow these mean 
field clusters to Bose condense. Actually it is well known 
among the nuclear physics community that even for such 
a small nucleus as the a-particle a momentum projected 
mean field approach yields a very reasonable description 
(TH - The reason for this stems, as already mentioned, 
from the presence of four different fermions, all attract- 
ing one another with about the same force. 

In Fig. [1] we sketch the situation, indicating that the 
two protons and two neutrons occupy the lowest OS" level 
of the mean field potential. Actually calculations show 
that the OS orbital of the self-consistent mean field resem- 
bles very much an oscillator wave function of Gaussian 
shape. In this respect the cartoon in Fig. [1] is not so 
far from reality. We suspect that the situation is generic 
for all strongly bound quartets which may be produced 
in the future and, therefore, our present study is of quite 
general interest. We will adopt this momentum projected 
mean field procedure in this work. 



II. THE IN-MEDIUM FOUR BODY EQUATION 

In-medium four body e quat ions are well documented 
in the literature since long [Ty] . In the present case of an 



in- medium quartet, the corresponding equation reads as 
follows [8|: 

(E — £l234) v I'l234 = (1 ~ /l ~ h) ^ V 12,V2' 1'2'34 

1'2' 

+ (1 - /l - /3) ^«13,1'3'*1'23'4 
1'3' 

+ permutations, (1) 

where ft = f( Si ) = [ e ( £ >-^/ T + l]" 1 with e. t = e{h) = 
kf/(2m) is the Fermi-Dirac distribution and £1234 = £1 + 
£2 + £3 + £4 (h = c = ks = 1: natural units). The matrix 
element of the interaction is i>i2,i'2' with the numbers 1, 
2, 3, • • • standing for all quantum numbers as momenta, 
spin, isospin, etc., as also in all other quantities in ([T]). 

In Eq. {!]), when E — A/i, this signals quartet conden- 
sation in very much the same manner as in the two body 
equation 

{E-e 12 )^i2 = (1 -fx - / 2 )^ t ' 12 4'2'*i'2', (2) 

1'2' 

where £12 = £1 + £2, the approach of T — > T c such that 
E — > 2/i signals the transition to a superconducting or 
superfluid state (the well known Thouless criterion [171]). 

Of course, as already stated several times, the determi- 
nation of T" needs the heavy solution of the in-medium 
modified four particle equation {T]). 

Following the discussion in the introduction, we, there- 
fore, make the following 'projected' mean field ansatz for 
the quartet wave function [J H [l8| , 

4 

*1234 = (27r) 3 ^ 3 )(k 1 +k 2 +k 3 + k 4 )[]( /3 (k l )x ST , (3) 

i=l 

where x is the spin-isospin function which we suppose 
to be the one of a scalar (S = T = 0). We will not 
further mention it from now on. We work in momen- 
tum space and <p(k) is the as-yet unknown single par- 
ticle OS wave function. In position space, this leads to 
the usual formula [HJ ^1234 — > / d 3 R J]^ =1 ^(r* - R) 
where fifa) is the Fourier transform of y(kj). If we 
take for (p(ki) Gaussian shape, this gives: I1234 — > 
cxp[— cJ2i<i<k<4:( r i ~ r fc) 2 ] which is the translationally 
invariant ansatz often used to describe a-clusters in nu- 
clei. For instance, it is also employed in the a-particle 
condensate wave function of Tohsaki, Horiuchi, Schuck, 
Ropke (THSR) in 0. 

Inserting the ansatz © into (Q]) and integrating over 
superfluous variables, or minimizing the energy, we arrive 
at the following non-linear, Hartree-Fock type of equation 
for the single particle OS wave function (p(k) = y(|k|) 

A(k)<p(k) + 3B(k) + 3C{k)(p{k) = 0, (4) 
where A(k), B(k), and C(k) are given by: 
r 4 w 3 t- T 4 P 

**->-/n££ Elk- 4 " 

J i=2 v ' Li=l 
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x |^ 2 )| 2 |^(fc 3 )| 2 |^(fc4)| 2 (27r) 3 ( 5( 3 )(^k 4 ) (5) 



1=1 



n000< i -/(«)-/(«» 



X «k 1 k 2 ,k' 1 k^(fci)¥'(fc2)¥'( fc 2)|^(A:3)| 2 |V5(fc4)| 2 
4 

x (2 7 r) 3 < 5( 3 )(^k l ) (6) 



»=i 



C(fci) 



n 000<>- /(«)-/(■»» 

x v k2k3M ' k/ w(k 2 )'p(k2)<f(k 3 )ip(k' 3 )\(p{k 4 )\ 2 



(7) 



where, in this pilot study, we neglect mean field shifts 
and effective mass contributions. 

From Eq. (f4]), we obtain the single particle wave func- 
tion in momentum space as 



<p(k) = 



-3B(k) 



A(k) + 3C(k) 



(8) 



As seen in Eqs. dU), {6]), and ([7|), since A(k), B(k), and 
C(k) depend on the wave function of ip(k), Eq. flSJ is 
strongly non-linear. Its solution can be found by itera- 
tion. 

For a general two body force Wk 1 k 2 ,k'k^) the equation 
to be solved is still rather complicated. We, therefore, 
proceed to the last simplification and replace the two 
body force by a unique separable one, that is 



Wkjka.k'jk^ 



where k = (ki - k 2 )/2, k' = (kj - k 2 )/2, K = ki + k 2 , 
and K' = k^ + k 2 . This means that we take a spin- 
isospin averaged two body interaction and disregard that 
in principle the force may be somewhat different in the 
S,T = 0, 1 or 1,0 channels. 

We are now ready to study the solution of {T]) for the 
critical temperature T C Q when the eigenvalue hits 4/x. For 
later comparison, the deuteron (pair) wave function at 
the critical temperature is also represented from Eqs. @ 
and |9j to be 



k 2 /m — 2/i 



d 3 k' 



k /k Q 



4>(k') 



( 10 ) 

where 4>{k) is a relative wave function of two parti- 
cles given by * 12 -> 0(| kl ~ k2 |)£ (3) (ki + k 2 ), and e = 
k 2 /(2m). From Eq. (fTtI|) . the critical temperature of pah- 
condensation is obtained with the following equation: 



' (27r) 3 fc 2 /m-2^ 



(11) 
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FIG. 2: Critical temperature of alpha and deuteron conden- 
sations as functions of chemical potential (a) and density (b), 
derived from Eq. ((4]) for the a-particle and Eq. for the 
deuteron, respectively. 



III. RESULTS FOR THE CRITICAL 
TEMPERATURE T c a 

In order to determine the critical temperature for a- 
condensation as a function of density n, we need to de- 
termine the chemical potential /z via 



n = 4 



d 3 k 



(12) 



and adjust the temperature so that the eigenvalue of ([I]) 
hits 4/j,. The two open constants A and fco in Eq- ® 
are determined so that binding energy (—28.3 MeV) and 
radius (1.71 fm) of the free ( /i = 0) a-particle come out 
right. The adjusted values are: A = —992 MeV fm 3 , and 
b = 1.43 fm -1 . The results of the calculation are shown 
in Fig. H 

In Fig. [21 the maximum of critical temperature T" max 
is at n = 5.5 MeV, and the a-condensation can exist up to 
Mmax = 11 MeV. It is very remarkable that the obtained 
results for T° well agree with a direct solution of (1) 
[l2| . These results for T° are by about 25 percent higher 
than the ones of our earlier publication We, however, 
checked that the underlying radius of the a-particle in 
that work is larger than the experimental value and that 
T" decreases with increasing radius of a-particle. Fur- 
thermore a different variational wave function was used 
in i. 

In Fig. [2] we also show the critical temperature for 
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deuteron condensation derived from Eq. (jlip . In this 
case, we take A = -1305 MeV fm 3 and k a = 1.46 fm -1 
to get experimental energy (—2.2 MeV) and radius (1.95 
fm) of the deuteron. It is seen that at higher densities 
deuteron condensation wins over the one of a-particlc. 
The latter breaks down rather abruptly at a critical pos- 
itive value of the chemical potential. Roughly speaking, 
this corresponds to the point where the a-particles start 
to overlap. This behavior stems from the fact that Fermi- 
Dirac distributions in the four body case, see ((T|), can 
never become step-like, as in the two body case, even not 
at zero temperature, since the pairs in an a-particle are 
always in motion. As a consequence, a-condensation gen- 
erally only exists as a BEC phase and the weak coupling 
regime is absent. 

Fig. [3] shows the normalized self-consistent solution 
of the wave function in momentum space derived from 
Eq. ijSJ) and the wave function in position space defined 
by its Fourier transform (p[r). Fig. [3]-(al) and (bl) are 
the wave functions of the free a-particle. As discussed in 
the Introduction, the wave function resembles a Gaussian 
and this shape is approximately maintained as long as fj, 
is negative, see Fig. [3]-(a2). On the contrary, the wave 
function of Fig. [5]-(a3), where the chemical potential is 
positive, has a dip around k — which is due to the Pauli 
blocking effect. For the even larger positive chemical po- 
tential of Fig. [5]-(a4) the wave function develops a node. 
This is because of the structure of the wave function de- 
rived in Eq. Q from where one can realize that again this 
stems from the Pauli blocking factor. The maximum of 
the wave function shifts to higher momenta and follows 
the increase of the Fermi momentum hp, as indicated 
on Fig. [3] From a certain point on the denominator in 
([8]) develops a zero and no self-consistent solution can be 
found any longer. 

On the other hand, the wave functions in position space 
in Figs. E2-(b2), (b3) and (b4) develop an oscillatory be- 
havior, as the chemical potential increases. This is remi- 
niscent to what happens in BCS theory for the pair wave 
function in position space p^ |. 



IV. DISCUSSION AND CONCLUSIONS 
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FIG. 3: Single particle wave functions in momentum space 
f>(k) (a), and in position space rtp(r) (b) at critical temper- 
ature, Eq. ((HJ. From top to bottom: (1) fi = -7.08 MeV, 
T c = MeV, n = fm" 3 (2) fx = -2.22 MeV, T c = 6.61 
MeV, n = 9.41 x 10~ 3 fin" 3 , (3) fj, = 6.17 MeV, T c = 8.45 
MeV, n = 3.07 x 1(T 2 fm 



3 , and (4) fx = 10.6 MeV, T c = 5.54 

3 



MeV, n = 3.34 x 10" 2 fm -3 . Figs, (al) and (bl) correspond 
to the wave functions for free Q-particle. The vertical lines in 
(a3) and (a4) are at the Fermi wave length &f = \/2mp,. 



In this work we again took up the study of the criti- 
cal temperature of a-particle (quartet) condensation in 
homogeneous symmetric nuclear matter. We essentially 
confirm the behavior of two previous studies [1, [l2j . The 
objective of the paper was to show that practically same 
results as before can be obtained with a strongly simpli- 
fying ansatz for the four particle wave function. Namely, 
this time, we used a momentum projected mean field 
variational wave function. This is based on the fact that 
the four different fermions of the quartet can occupy the 
same single particle OS'-wave function in the mean field. 
The latter is to be determined from a self-consistent non 
linear HF-type of equation as a function of chemical po- 
tential or density. The relation between the chemical 



potential and density is taken from the free Fermi gas 
relation, Eq. (| 1 2|) . However, the total nucleon density 
of the system must be calculated from the mean single 
nucleon state occupation number taking into account cor- 
relations, so that the contribution of bound states to the 
total nucleon density is taken into account, see Ref. [20] |. 
To calculate the critical temperature not as function of 
the free nucleon density, see Fig. [2]b, but of the total 
nucleon density, a generalization a la NSR [l3[ must be 
performed, that is we have at least to incorporate the 
contribution of the a-particle density including the con- 
densate to the single particle occupation numbers. This 
shall be investigated in future work. 
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Besides, in this work, we used the isospin-independent 
separable potential, Eq. © , for the two-body interaction 
as a simplification. Comparison with a realistic two-body 
interaction is interesting. This study also shall be done 
in the future. 

The self-consistent wave function has been studied in 
momentum and position space. For negative chemical 
potential the single particle wave function behaves like 
a Gaussian. However, once the chemical potential turns 
positive, then the single particle wave function in r-space 
starts to oscillate. This is a well known feature from 
ordinary pairing. 

We, therefore, have demonstrated that a very simpli- 
fying momentum projected mean field ansatz suffices to 
account for the salient features of quartet condensation. 
This is very helpful for the next step which is more com- 
plicated, i.e. the incorporation of quartet condensation 
self-consistently into the Equation of State (EOS). 

We should, however, be aware of the fact that our pro- 
jected mean field ansatz for the quartet wave function 
can only be a valid approximation so long as well de- 
fined quartets exist. In the break down region seen on 
Fig. [21 this is certainly no longer the case. How the quar- 
tet phase evolves into a superfluid phase of pairs is an 
open question. A possibility to study this very inter- 
esting problem could be to write down the in-medium 
four body equation JT]) directly in the BCS formalism, 
i.e. with the corresponding BCS coherence factors. It 
may be foreseen that the latter only catch on close to 



the transition region. Another interesting problem for 
the future is how the present results are modified in the 
asymmetric case, that is in the case of neutron excess. 

The success of our study to employ a very simplifying 
ansatz of the mean field type for the quartet wave func- 
tion, may open wide perspectives. Besides to push the 
description of quartet condensation much further, there 
might exist the possibility that even for the case of a gas 
of trions such a projected mean field ansatz is a quite 
valid approach. In the case of three colors, like quarks 
in the constituent quark model for nucleons, a harmonic 
confining potential is frequently assumed and the three 
quarks can occupy the lowest 05* state, analogously to the 
case of quartets treated in the present paper. Of course, 
trions are composite fermions and cannot be treated in 
the same way as bosonic composites, since they shall form 
a new Fermi gas with their own new Fermi level. How 
this situation can eventually be treated has recently been 
outlined in [2l| . 
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